Double symmetry breaking of solitons in one-dimensional virtual photonic crystals 
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We demonstrate that spatial solitons undergo two consecutive spontaneous symmetry breakings 
(SSBs), with the increase of the total power, in nonlinear photonic crystals (PhCs) built as arrays 
of alternating linear and nonlinear stripes, in the case when maxima of the effective refractive index 
coincide with minima of the self-focusing coefficient, and vice versa, i.e., the corresponding linear and 
nonlinear periodic potentials are in competition. This setting may be induced, as a virtual PhC, 
, by means of the EIT (electromagnetically-induced-transparency) technique, in a uniform optical 

■ medium. It may also be realized as a Bose-Einstein condensate (BEC) subject to the action of 

combined periodic optical potential and periodically modulated Feshbach resonance. The first SSB 
happens at the center of a linear stripe, pushing a broad low-power soliton into an adjacent nonlinear 
^ft, . stripe and gradually suppressing side peaks in the soliton's shape. Then, the soliton restores its 

' symmetry, being pinned to the midpoint of the nonlinear stripe. The second SSB occurs at higher 

' powers, pushing the narrow soliton off the center of the nonlinear channel, while the soliton keeps its 

£N| . internal symmetry. The results are obtained by means of numerical and analytical methods. They 

may be employed to control switching of light beams by means of the varying power. 
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I. INTRODUCTION 

It is commonly known that linear eigenstates supported by symmetric potentials, in contexts such as quantum 
\Jq , mechanics and photonic crystals (PhCs), may be classified according to representations of the underlying symmetry 
£ — , ■ group [l| . The addition of the nonlinearity frequently gives rise to the effect of the spontaneous symmetry breaking 
t-H \ (SSB), i.e., reduction of the full symmetry group to its subgroups, in the generic situation. The basic feature of 
. the SSB is the transition from the symmetric ground state (GS) to that which does not follow the symmetry of the 
' potential. The simplest manifestations of the SSB, which were predicted in early works Q, and then in the model 
of nonlinear dual-core optical fibers 13, 0], occur in settings based on symmetric double- well potentials, or similarly 
t-H ' designed nonlinear pseudopotentials [5| (the latter name for effective potentials induced by a spatially inhomogeneous 
nonlinearity is used in solid-state physics @). In the quantum theory, the nonlinearity naturally appears in the 
context of the mean- field description of Bose-Einstein condensates (BECs), with the Schrodinger equation replaced 
by the Gross-Pitaevskii equation [?J • Similarly, the self- focusing makes PhCs a medium in which the linear symmetry 
competes with the nonlinearity, in two-dimensional (2D) Q and ID settings. 

A natural extension of the SSB in the double- well potential is the transition from symmetric to asymmetric solitons 
in the geometry which adds a uniform transverse dimension to the potential, extending the potential from the double 
well into a double trough. In this setting, the soliton self-traps in the free direction due to the self-focusing nonlinearity. 
The SSB effect in the solitons may be described by means of the effectively ID two-mode approximation, which is 
tantamount to the usual temporal-domain model of dual-core o ptic al fibers (3] , and may also be applied to the BEC 
loaded into a pair of parallel tunnel-coupled cigar-shaped traps [13|. In the application to the BEC, a more accurate 
description of the symmetry-broken solitons was developed too, in the framework of the 2D Gross-Pitaevskii equation, 
for the linear [l4[ and nonlinear [l5| double-trough potentials. 

PhC media are modeled by combinations of linear and nonlinear potentials, which correspond to the alternation 
of material elements and voids in the PhC structure A similar setting, emulating the PhC,may be induced 

by means of the clcctromagnetically-induced-transparency (EIT) technique in a uniform medium [lfjj . For BEC, a 
counterpart of the PhC may be generated as a combination of the linear potential, induced by the optical lattice, 
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FIG. 1: (Color online) (a) The scheme of the virtual photonic crystal with period D — rfi + d,2- The blue (darker) and 
gray (lighter) slabs represent the nonlinear and linear stripes, respectively, (b) The corresponding Kronig-Penney modulation 
function. 



and a nonlinear pseudopotential imposed by a periodically patterned magnetic or optical field modifying the local 
nonlinearity via the Feshbach-rcsonancc effect. Actually, the latter setting may be more general than the PhC, as the 
so designed linear and nonlinear potentials may be created with incommensurate periodicities fl7j . 

Solitons in periodic linear and nonlinear potentials have been studied theoretically in many works, as reviewed 
in Refs. [l8| and [lj|. In particular, specific types of gap solitons were predicted in ID models of PhCs featuring 
the competition between the periodically modulated refractive index and self-defocusing material nonlinearity [l2| . 
Spatial optical solitons, supported by effective linear potentials, were created in various experimental setups [201 ] - 

Unlike the double-well settings, periodic potentials usually do not give rise to SSB in solitons, although examples 
of asymmetric solitons were found in some ID models (Tlj . Indeed, ID optical media built as periodic alternations of 
self-focusing material elements and voids feature no competition between the effective linear and nonlinear potentials, 
as minima and maxima of both types of the potentials coincide, hence there is no drive for the SSB in the medium 
(as mentioned above, the competition takes place if the material nonlinearity is self-defocusing, but in that case the 
corresponding gap solitons do not feature SSB either [l2| )■ Competing potentials leading to SSB effects might be 
possible if maxima of the refractive index would correspond to minima of the self-focusing nonlinearity. While this is 
impossible in usual PhCs, composed of material stripes separated by voids, the effective potential structures induced 
by EIT patterns in uniform media admit such a situation [l6j (as said above, a similar setting is also possible in 
BEC [13 EH)- The objective of this work is to study the SSB for spatial solitons in the virtual PhC of this type, 
following the gradual increase of the total power of the soliton. It will be demonstrated that the symmetry breaking 
happens twice in this setting, first to a low-power soliton centered around a midpoint of a linear channel, and then to 
a high-power beam situated at the center of the nonlinear stripe. 

The setting under the consideration is displayed in Fig. HJa), obeying the following version of the nonlinear 
Schrodinger equation for local amplitude V(x, z), which is a function of propagation distance z and transverse coor- 
dinate x: 



iu z = -(l/2)u xx + V(x) (1- \u\ 2 )u, (1) 

where the Kronig-Penney (KP) potential function, V(x), is defined as per Fig. [ljb). We stress that the self-focusing 
sign of the nonlinearity makes the model different from ones with the competition between the linear and nonlinear 
potentials provided by the self-defocusing [l2j. We consider the version of the system with d\ = c?2 = d in Fig. Ufa). 

The scaling is fixed by setting the modulation depth to be Vo = 0.02, which leaves half-period d of the potential as 
a free parameter. In Section 2, we present numerical results, at first, for d — 10, which illustrates a generic situation. 
Then, we demonstrate the results for other values of the period, namely, d — 5,8, 11, and 14. In particular, it will be 
demonstrated that the second SSB vanishes at large values of d. In Section 3, we report analytical approximations, 



II. SOLITONS AND THEIR SYMMETRY BREAKINGS (NUMERICAL RESULTS) 



A. Results for d — 10 



GS (ground-state) solutions to Eq. (p} with the above-mentioned values of the parameters, d\ = c?2 = 10 and 
Vo = 0.02, were found by means of the imaginary-time-propagation method (2lj . Characteristic features of GS solitons, 
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FIG. 2: (Color online) Profiles of ground-state solitons found at the following values of the total power: P = 12 (a), P = 16 
(b), P = 29 (c), and P = 50 (d). 



which distinguish them from bound complexes of two or several solitary beams, is the presence of a pronounced central 
peak, and the absence of nodes (zero crossings). 

Figure [2] displays a characteristic set of GS profiles found at different values of total power, which is defined as 
P = f-°° \u(x)\ 2 dx. Representing the stationary solutions as u(x,z) = e~ tfJ,z U(x), it is straightforward to check 
that wavenumber — fj, of all the GS solitons falls into the semi-infinite gap, in terms of the spectrum generated by 
the linearized version of Eq. (JTJ) , see Fig. |3Jd) below; this situation is natural for the system with the self- focusing 
nonlincarity. Further, real-time simulations of Eq. (fTJ) (not shown here) demonstrates that all the GS modes are 
stable against perturbations. As seen from Fig. Eld), the stability of the GS solitons is also supported by the 
Vakhitov-Kolokolov criterion [23j], d/i/dP < 0. 

At sufficiently small values of P (in the weakly nonlinear regime), Fig. HJa) demonstrates that the central and side 
peaks of the GS soliton are situated in linear channels, the soliton being symmetric about its central peak, which 
exactly coincides with the midpoint of the corresponding linear stripe. Although the light is chiefly guided by the 
linear stripes in this regime, the soliton of course cannot exist without the self-focusing nonlincarity, even if it is weak. 

As seen in Fig. H^b), the first SSB event happens with the increase of P, breaking the symmetry of the weakly- 
nonlinear GS soliton and spontaneously shifting its central peak off the center of the linear channel in which the peak 
is trapped. The symmetry of the side peaks gets broken too, although they remain trapped in the linear channels. It 
is relevant to stress that as ym metric GS solitons were not reported in previously studied versions of nonlinear systems 
with KP potentials fill. IT2I. 22 1, although some exact asymmetric solutions for non- fundamental solitons, which feature 
nodes, were found in Ref. [11 1. 

With the further increase of the power, the GS soliton undergoes strong self-compression, which eliminates all side 
peaks, while the central peak moves from the linear stripe into an adjacent nonlinear one, ending up in the middle of 
the nonlinear stripe. A seen in Fig. HJc), in the corresponding moderately nonlinear regime the single-peak soliton 
eventually restores its symmetry about the midpoint of the nonlinear channel into which it has shifted. Accordingly, 
light is guided by the nonlinear stripe in this regime. 

Finally, in the strongly nonlinear regime the second SSB event happens, as seen in Fig. [U[d), where the narrow 
GS soliton spontaneously shifts from the midpoint of the nonlinear stripe, although staying in it. To the best of our 
knowledge, such a repeated SSB of solitons has not been reported in other models of nonlinear optics or BEC. The 



4 




10 15 20 25 30 35 40 45 50 10 15 20 25 30 35 40 45 50 

P P 



(c) (d) 

FIG. 3: (Color online) (a) The coordinate of the soliton's center of mass versus total power P. (b) The duty cycle (DC), showing 
the share of the power trapped in linear stripes, and the soliton's asymmetry measure (AS), versus P. (c,d): The soliton's width 
(W a ) and propagation constant (— /i) as functions of P (it can be checked that all values of /i belong to the semi-infinite gap, 
in terms of the spectrum of the linearized system). Note that the n(P) dependence satisfies the Vakhitov-Kolokolov stability 
criterion [H, d^/dP < 0. 



second SSB seems a counter-intuitive effect, as one might "naively" expect that, with indefinite increase of total power 
P, the narrow high-power soliton would be only stronger nested at the center of the nonlinear stripe. Nevertheless, 
an explanation to this effect is possible, as argued below. 

To quantify the double SSB, we define the following characteristics of the soliton, as functions of the total power, 
P: center-of-mass coordinate x mc , average width W a , linear-stripe duty cycle (DC), the soliton's asymmetry measure 
(AS), and the above-mentioned propagation constant, —\x (if the system is realized as the BEC model, // is the 
chemical potential): 



/+oo 
x\u\ 2 dx, 
-oo 



DC = P- 1 / \u\ 2 dx, AS = P~ 

Jhs Jo 



[u(x max - y)] 2 - [u(x max + y)\ 



f+oo /> + oo 

r2 d— 1 / /„ ™ \2i ,|2 j ,, _ d— 1 



dy, 

w l = P ~ X I (x~ x mc ) 2 \u\ 2 dx, h^P- 1 I u*Hudx, (2) 

J — oo J — oo 

where J Ls stands for the integral taken over the linear stripes, x max is the location of the maximum value of |u(x)|, 
and the Hamiltonian operator is H = (l/2)d xx + V(x)[l — \u\ 2 ]. For DC > 50% (< 50%), light is mainly guided by 
the linear (nonlinear) stripes. The strength of the SSB as a whole is quantified by shift x mc , while AS quantifies the 
related inner symmetry breaking of the soliton. 

The overall description of the double SSB is provided, in Fig. [3l by plots showing the evolution of quantities ([2]) 
with the increase of P. In panel (a), the dashed lines, x = and x = —10, mark the positions of the two symmetric 
points in the KP potential, which correspond, respectively, to midpoints of the linear and nonlinear stripes. The plots 
in Figs. |3Ja,b) clearly demonstrate that the low-power GS soliton remains symmetric, being centered at x — 0, for 
P < 15. 
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FIG. 4: (Color online) x ulc (P) for different value of the lattice period: (a) d — 5, (b) d = 8, (c) d = 11, and (d) d = 14 



The first SSB event occurs at Pggg ~ 15. The soliton becomes asymmetric, gradually moving from the midpoint of 
the linear stripe (x = 0) towards the center of the adjacent nonlinear one (x — —10), in the interval of 15 < P < 22. 
As seen in Fig. [3fb) , the soliton attains the largest inner asymmetry degree, AS > 40%, around P ~ 18. Observe that 
the asymmetric soliton gradually loses its side peaks, simplifying into the single-peak shape, as seen in Figs, dta-c). 

The dependence AS(P) in Fig. GUb) shows that the SSB occurring at P = Pggg is of the supercritical type, i.e., it 
may be identified with a phase transition of the second kind 24 1 . 

Next, the moderately high-power soliton remains completely symmetric, centered at the midpoint of the nonlinear 
channel, at 22 < P < 30. Note also that the duty cycle falls to values DC < 50% at P > 22, which implies the switch 
from the quasi-linear guidance to that dominated by the nonlinear pscudopotcntial. 

(2) 

The second SSB occurs at P^ S ' B sa 30. At P > 30, the high -power soliton gradually shifts from x = — 10 towards 
the edge of the nonlinear stripe, while keeping a virtually undisturbed symmetric shape, with AS = 0, see Fig. [3]Jb). 
This instance of the SSB may also be considered as a phase transition of the second kind. 



B. Extension to other values of the potential's period 



To present the most general results, in Fig. 0] we report dependence x 
d\ = c?2 = d), viz., d = 5, 8, 11, and 14. We notice that, when d increases from 5 in Fig 



(2) 

power at which the second SSB takes place decreases from Pggg ~ 120 to ~ 24. This is explained by the fact that 
the second SSB requires the average width of the soliton to be much narrower than the width of the nonlinear stripe, 
hence, with smaller d, larger P is needed to sustain the second SSB for a tighter localized soliton. This argument 
also explains the shrinkage of the interval in which the soliton is pinned to the symmetric position at the midpoint 
of the the nonlinear stripe with the increase of d, as seen in Figs. 4(a)| (b) and (c). If d keeps increasing, this interval 
eventually disappears at d > 11, as shown in Fig. |4(d)| pertaining to d — 14. 

Further, in Fig. [S]we plot soliton characteristics DC, AS and W a for d = 14. It exhibits a direct transition from 
the original symmetric state to the final asymmetric state, without the second SSB. 



4(a) to 11 in Fig. 



(P) for different values of d (still with 

4(g) the 



6 





(a) 



(b) 



FIG. 5: (Color online) (a) Dependences DC(P) and AS(P) for d = 14. (b) VK a (P) for d = 14. 



III. ANALYTICAL CONSIDERATIONS 

Both instances of the SSB revealed by the numerical findings can be explained by means of analytical approxima- 
tions. To address the first SSB, which happens to the broad low-power GS soliton, we replace the KP modulation 
function with period D = 2d in Eq. ([T]). V(x) (see Fig. [!}, by the combination of its mean value and the first harmonic 
component, dropping higher-order harmonic components (cf. Ref. (Toj): 



Accordingly, Eq. (JTJ is replaced by 



V(x) = (V /2) [1 - (4/tt) cos (2ttx/D)}, 



iu z = -(l/2)u xx + V(x) (1 - M 2 ) u. 



(3) 



(4) 



In the zero-order approximation, one may neglect the variable part of the modulation function, whose period is much 
smaller than the width of soliton. In this case, Eq. ^ gives rise to the obvious soliton solution, 



u (x, z) = ^/27VbV fc *sech {r, (x - £)) , k = \ (rr 2 - V Q ) , 

with coordinate of the soliton's center £ and inverse width 77. The total power of the soliton is P = Arj/Vo. 
The Hamiltonian corresponding to Eq. ((4]) is 



(5) 



H = - 

2 



-j-oc 



V(x) [2| 



dx. 



(0) 



As follows from here, the energy of the interaction of the broad soliton, taken as per Eq. (JSJ) (which neglects the 
distortion of the soliton's shape under the action of the potential) with the variable part of modulation function ([3]) , 
expressed in terms of the soliton's power, is 



U(0 = - (VoA) 



cos(2ttx/£>) 2\u(x)\ 2 - \u(x)f 



dx 



8tt 
~~d 



sinh 



/ 4tt 2 



\PV D 



-1 r 



-1 



2 
3^ 



7T \ 2 
Dj 



16 



(7) 



For small P (low-power solitons), energy ([7]) gives rise to a local minimum of the corresponding potential, i.e., a stable 
equilibrium position of the soliton at £ = 0, provided that 



V Q > 2ir 2 / (3D 2 



(8) 



Note that the parameter values adopted above to generate Figs. [TJ21 D = 20 and Vo = 0.02, satisfy this condition. 
Then, with the increase of P, the equilibrium position predicted by potential ([7]) becomes unstable at 



P > P< 



jg ee (4/V ) ^(3/2)Vo-(tt/D) 2 



(9) 
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At P > P$ B , the soliton moves away from £ 
and D = 20, expression (J9j> yields Pgg B 



above from the numerical data, P 



(i) 



0, i.e., this is the point of the first SSB. Note that, for Vo = 0.02 

i 14.6, which practically exactly coincides with the first SSB point identified 

15. If D is too small and does not satisfy condition © for given Vo, this 
means that the simplest approximation cannot be used, and corrections to the average form of the soliton ([5]) should 
be taken into account, which is beyond the scope of the present analysis. 

Proceeding to the analytical consideration of the second SSB, we rewrite Eq. (JTJ as 



SSB 



iu z = -(l/2)u m -V(y)(\u\ 2 -l)- 



(10) 



where coordinate y is defined so as to place the center of the nonlinear stripe at y = 0. Aiming to consider narrow 
solitons, trapped in the given channel, which do not feel the presence of other nonlinear stripes, we define the 
modulation function here so that V{y) = Vq in the nonlinear stripe, and V{y) = outside of it. Then, the narrow 
soliton with the center located at point y = £ (generally, shifted off the center of the nonlinear stripe) has the following 
form: 



u (z, x) 



ikz 



! exp 



'd% 
* 2 



sech (ji (y 
£ sgn (y)) - 



0), at \y \ < 
VV 2 - 2F (M 



<1± 

2 



at \y\ > 



(11) 



where this time the inverse width of the soliton is assumed to be large, r\ 3> l/cfe, k is the same as in Eq. ([5]), and 
the total power of the narrow soliton is P ps 2?y/Vb. 

The substitution of the wave field (fTTj) into the Hamiltonian of Eq. (|TU)) yields the following effective potential of 
the interaction of the soliton with the nonlinear stripe which holds it: 



U(£) = -2ri 




e- ,,d2 cosh (2r?0 + l^e" 2 '^ 2 cosh (4t/£) 
Vo 



(12) 



cf. Eq. 0. As might be expected, the last term in this potential tends to keep the soliton at the center (£ = 0), while 
the other terms push it towards the edge of the nonlinear stripe. The equilibrium position is defined by equation 
dll/dt; = 0. The substitution of potential (jT2"j) into this equation yields two solutions: either £ = 0, which corresponds 
to the soliton placed exactly at the center, and the off-center equilibrium, determined by the following expression: 



cosh (2?7£) = i 

Solution (|13|) exists if it yields cosh (2?7£) > 1, i.e., 



o.Vd-2 



(13) 



" d2 >8 1- ^l-2Vo/v 2 ) +2V /r) 



(14) 

Following the assumption that the soliton is narrow, we assume 2Vo/?? 2 1, hence Eqs. (fT3"|) and (H3J) are simplified 
as follows: 



cosh (V P0 - 3 {2V P 2 ) 1 exp (Vbd 2 P/2) . 
exp (VodaP/2) > (2/3)V P 2 . 



(15) 
(16) 



The second SSB is realized, in the framework of the present approximation, as the displacement of the soliton 
from £ = to point ([TB")) . hence inequality (fTo| . if replaced by the respective equality, offers a rough approximation 

(2) 

for the second SSB point, Pg SB ["rough" because it was derived taking into regard exponentially small terms in Eq. 
(fT2"|) . which is a crude but meaningful approximation [25|]. In particular, Eq. (fTB")) predicts that, for Vo = 0.02, the 

(2) 

numerically found value reported above, Pg SB = 30, corresponds to cfe — 8.3, which is not far from d2 = 10 which was 
actually used. For P — > oo, Eq. (|15p yields £ — > ±^2/2, i.e., within the framework of the present approximation, the 
soliton moves to the edge of the nonlinear stripe. On the other hand, for larger values of c?2 inequality (|16p always 
holds, which explains the disappearance of the second SSB in the numerical picture displayed in Fig. |4j 



IV. CONCLUSION 



We have demonstrated the effect of the double SSB (spontaneous symmetry breaking) for GS (ground-state) stable 
solitons in the ID medium with the competing periodic linear potential and its nonlinear counterpart (pscudopotcntial) 
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induced by a periodic modulation of the local self-attraction coefficient. This medium may be realized as a virtual 
PhC (photonic crystal) imprinted by means of the EIT technique into a uniform optical medium, and also as the 
BEC setting using a combination of an optical lattice and the spatially periodic modulation of the nonlinearity via 
the Feshbach resonance. The two SSB events occur in the low- and high-power regimes, pushing the soliton off the 
symmetric positions at the center of the linear and nonlinear stripes, respectively. In the former case, the SSB also 
affects the shape of the low-power soliton, making it asymmetric and gradually stripping it of side peaks. In the latter 
case, the narrow high-power soliton, while shifting off the midpoint of the nonlinear channel, keeps the symmetric 
shape. At intermediate values of the power, the soliton is completely symmetric, staying pinned at the center of the 
nonlinear stripe. On the other hand, the increase of the period of the potential structure leads to the direct transition 
from the original symmetric state to the final asymmetric one, while the second SSB point disappears. These results 
, which were obtained by means of systematic numerical computations and explained with the help of analytical 
approximations, suggest a possibility to control the switching of spatial optical solitons in the virtual PhC by varying 
their power. 

This work may be extended in other directions. It particular, a challenging problem is to investigate similar settings 
and effects for 2D solitons, in terms of PhCs and BEC alike. 
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by the Chinese agencies NKBRSF (grant No. G2010CB923204) and CNNSF(grant No. 10934011). 



[1] L. D. Landau and E. M. Lifshitz, Quantum Mechanics (Moscow: Nauka Publishers, 1974). 

[2] E. B. Davies, Comm. Math. Phys. 64, 191 (1979); J. C. Eilbeck, P. S. Lomdahl, and A. C. Scott, Physica D 16, 318 (1985). 

[3] A. W. Snyder, D. J. Mitchell, L. Poladian, D. R. Rowland, and Y. Chen, J. Opt. Soc. Am. B 8, 2101 (1991). 

[4] C. Pare and M. Florjahczyk, Phys. Rev. A 41, 6287 (1990); A. I. Maimistov, Kvant. Elektron. 18, 758 [Sov. J. Quantum 
Electron. 21, 687 (1991)]; N. Akhmediev and A. Ankiewicz, Phys. Rev. Lett. 70, 2395 (1993); P. L. Chu, B. A. Malomed, 
and G. D. Peng, J. Opt. Soc. A B 10, 1379 (1993); B. A. Malomed, in: Progr. Optics 43, 71 (E. Wolf, editor: North 
Holland, Amsterdam, 2002). 

[5] L. C. Qian, M. L. Wall, S. Zhang, Z. Zhou, and H. Pu, Phys. Rev. A 77, 013611 (2008); T. Mayteevarunyoo, B. A. 

Malomed, and G. Dong, ibid. A 78, 053601 (2008). 
[6] W. A. Harrison, Pseudopotentials in the Theory of Metals (Benjamin: New York, 1966). 
[7] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Clarendon Press: Oxford, 2003). 

[8] P. Xie, Z.-Q. Zhang, and X. Zhang, Phys. Rev. E 67, 026607 (2003); A. Ferrando, M. Zacares, P. Fernandez de Cordoba, 
D. Binosi, and J. A. Monsoriu, Opt. Exp. 11, 452 (2003); 12, 817 (2004); J. R. Salgueiro, Y.i S. Kivshar, D. E. Pelinovsky, 
V. Simon, and H. Michinel, Stud. Appl. Math. 115, 157 (2005); A. S. Desyatnikov, N. Sagemerten, R. Fischer, B. Terhalle, 
D. Trager, D. N. Neshev, A. Dreischuh, C. Denz, W. Krolikowski, and Y. S. Kivshar, ibid. 14, 2851 (2006). 

[9] Q. Li, C. T. Chan, K. M. Ho and C. M. Soukoulis, Phys. Rev. B 53, 15577 (1996); E. Lidorikis, Q. Li, and C. M. Soukoulis, 
ibid. 54, 10249 (1996). 

[10] B. A. Malomed, Z. H. Wang, P. L. Chu, and G. D. Peng, J. Opt. Soc. Am. B 16, 1197 (1999). 
[11] Y. Kominis, Phys. Rev. E. 73, 066619 (2006); Y. Kominis and K. Hizanidis, Opt. Exp. 16, 12124 (2008). 
[12] Y. Kominis and K. Hizanidis, Opt. Lett. 31, 2888 (2006); T. Mayteevarunyoo and B. A. Malomed, J. Opt. Soc. Am. B 
25, 1854 (2008). 

[13] A. Gubeskys and B. A. Malomed, Phys. Rev. A 75, 063602 (2007); A 76, 043623 (2007). 

[14] M. Matuszewski, B. A. Malomed, and M. Trippenbach, Phys. Rev. A 75, 063621 (2007); M. Trippenbach, E. Infeld, J. 

Gocalek, M. Matuszewski, M. Oberthaler, and B. A. Malomed, ibid. A 78, 013603 (2008). 
[15] N. V. Hung, P. Ziri, M. Trippenbach, and B. A. Malomed, Phys. Rev. E 82, 046602 (2010). 

[16] C. Hang and V. V. Konotop, Phys. Rev. A 81, 053849 (2010); Y. Li, B. A. Malomed, M. Feng, and J. Zhou, ibid. 82, 
633813 (2010). 

[17] H. Sakaguchi and B. A. Malomed, Phys. Rev. A 81, 013624 (2010). 

[18] Y. V. Kartashov, V. A. Vysloukh, and L. Torner, in: Progr. Optics 52, 63 (E. Wolf, editor: North Holland, Amsterdam, 
2009). 

[19] Y. V. Kartashov, B. A. Malomed, and L. Torner, Solitons in nonlinear lattices, Rev. Mod. Phys., in press. 

[20] F. Lederer, G. I. Stegeman, D. N. Christodoulides, G. Assanto, M. Segev, and Y. Silberberg, Phys. Rep. 463, 1 (2008). 

[21] M. L. Chiofalo, S. Succi, and M. P. Tosi, Phys. Rev. E. 62, 7438 (2000). 

[22] I. M. Merhasin, B. V. Gisin, R. Driben, and B. A. Malomed, Phys. Rev. E. 71, 016613 (2005). 

[23] M. Vakhitov and A. Kolokolov, Radiophys. Quantum. Electron. 16, 783 (1973); L. Berge, Phys. Rep. 303, 259 (1998). 

[24] L. D. Landau and E. M. Lifshitz, Statistical Physics (Moscow: Nauka Publishers, 1976). 

[25] Yu. S. Kivshar and B. A. Malomed, Phys. Rev. Lett. 60, 164 (1988); Rev. Mod. Phys. 61, 763 (1989). 



